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Abstract 

In order to study the structure of excited states we perform a variational 
calculation after spin parity projection (VAP) within the framework of An- 
tisymmetrized Molecular Dynamics (AMD). The framework is proven to be 
a new powerful approach for the study of the various structures of excited 
states because it is free from model assumptions such as inert cores, existence 
of clusters, and the axial symmetry. By using finite range interactions with a 
density dependent term we reproduce well all the energy levels below 15 MeV 
in 12 C. This is the first theoretical model that reproduces many E2 transition 
rates and [3 decays to 12 C successfully. 
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Clustering is one of the important features of the light nuclei. The Ikeda diagram |J 
predicts molecule-like structures in the excited states, according to the activation of the 
inter-cluster relative motion. The cluster model has been successful for the study of the 
molecule-like states. However the model is not suitable for describing the shell-model-like 
features, such as the j-j coupling and the particle-hole excitations. The model, as we 
see later, failed to reproduce important data such as the level spacing and the transition 
strength because the dissociation of clusters is not taken into account by the conventional 
cluster model. These shell-model-like aspects are often found in many states including the 
ground state. On the other hand, the shell model successfully reproduces properties of 
many excited levels of light nuclei (for example but a number of states in light nuclei 

have been left unsolved because it is difficult for the shell model to describe well-developed 
molecule-like states which may exist in low energy region. 

There are very few approaches which can describe both the molecule-like and the shell- 
model-like phenomena. For instance, in 12 C which has been studied for a long time , 
the structure of the excited O^, O3 , and 2% states is not yet clarified. Therefore, it is desirable 
to develop a theoretical approach that is able to deal with both aspects in a unified manner. 
Our aim is to systematically study how the structure of nuclei changes with increasing 
excitation energy. 

The AMD has already proved to be a very useful theoretical approach for the structure 
of general light nuclei [3,R |l0l , |lll . The AMD basis wave functions are written as Slater 
determinants where the spatial part of each single particle wave function is a Gaussian. An 
important point of AMD is that we do not need model assumptions such as inert cores and 
the axial symmetry. In the previous works on the ground states of light unstable nuclei 



P JlQjTT|| , we applied a simplest version of AMD in which we make energy variation only 
after parity projection of a Slater determinant. The simplest version of AMD was found 
to describe the properties of ground states successfully. We also studied the change in 
structure along the yrast line of 20 Ne by a cranking method of AMD M instead of the 
angular momentum projection before energy variation. 

In this letter, we propose a new microscopic approach of variational calculation after 
parity and total angular momentum projection within the framework of antisymmetrized 
molecular dynamics (AMD). In order to confirm the userfulness of the method, we apply this 
approach for the first time to the study of excited states of 12 C using finite range interactions. 

Concerning the formulation of AMD, the reader is refered to papers ||,|9],[10| • An AMD 
wave function is a Slater determinant of Gaussian wave packets; 

®amd = ^A{ipx,p 2 ,- • ■ ,(p A }, (1) 



<Pi = faiXiTi '■ ' 




(2) 



where the centers of Gaussians Zj's are complex variational parameters. Xi is the intrin- 
sic spin function parametrized by and is the isospin function which is fixed to be 
up(proton) or down(neutron) in the present calculation. We vary the energy expectation 
values for the parity and total angular momentum eigenstates (YAP calculation), 
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where the operator of total angular momentum projection P( IK i is / dQD^ K ,(Q)R(Q). The 
expectation values are calculated numerically by the sum of mesh points instead of integral 
on Euler angle fl We make energy variation by the use of the frictional cooling method j9]. 
Here we represent the intrinsic wave function <&amd obtained by VAP for the lowest states 
with a given spin parity as . Higher excited states are constructed by superposing 
wave functions so as to orthogonalize to the lower states. That is to say that the n-th 
state is calculated by varing the energy expectation value of the othogonal component 
to the lower states; 

pfic*? - eWwi^w^- (4) 

fc=l 

By making VAP calculations with various sets of {J±,n} we obtain many intrinsic states 
{$!,-• -,$ m }, which approximately correspond to the states. Final results are con- 
structed by diagonalizing the hamiltonian matrix (Pj^,$j|i/|P/^,,$,,) formed from all the 
intrinsic states. The resonance states are treated within a bound state approximation. 

In this work the adopted interactions are the central force of the modified Volkov nol of 
case 3 |T2j which consists of the finite range two-body force and the zero range three-body 
repulsive force as a density dependent term, the spin-orbit force of G3RS with two range 
Gaussians [[HJ and the Coulomb force. The Majorana parameter used here is m = 0.62, and 
the strength of G3RS force is U\ = —u^ = 3000 MeV. We chose the width parameter v for 
Gaussians in Eq.|| as 0.18 . 

Energy levels of 12 C are shown in Figjl] compared with the experimental data. The 
AMD calculations reproduce the energy levels very well. In the cluster model calculations 
the level spacing between 0^ and 2\ was always underestimated. For instance, in a GCM 
calculation || it was 2.2 MeV which is much smaller than the experimental value 4.4MeV. 
In the present result the large level spacing between 0^ and 2\ agree with the experimental 
data. We think that it is because the theory describes successfully the dissociation of a 
cluster in the ground 0^ state due to the LS force. It is found that molecule-like states with 
well-developed 3a structures construct rotational bands of K n =3~, 1~, Of, Of. Even by a 
(0 + 2)hu) shell-model calculation for p — shell nuclei they can not describe Of and 2f 
states in the K n = Of band because such states contain highly excited n Tiu components in 
terms of the shell model. 

The transition strengths are of great help to investigate the structure of excited states. 
The theoretical and the experimental values of E2 transition strengths are shown in Table |. 
In a GCM calculation || of the 3a cluster model the B(E2; 0^ — * 2^) was about one third 
of the experimental data (see Table Q. The present result for B(E2; Of — > 2f) is as much as 
the experimental data. 

The calculated strengths of Gamov- Teller j3 decay transitions from 12 N to the excited 
states of 12 C are shown in Table |T| and compared with the experimental data. The ground 
state of the parent nucleus 12 N is obtained with a variational calculation after spin parity 
projection to 1 + states. It is difficult to discuss (3 transitions to 12 C* within the framework of 
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TABLE I. The strength of the E2 transitions of 12 C. 



Transition 


exp. 


present 


GCM | 


I 


B(E2; 2+ -+0+) 
B(E2; 0+ ^2+) 


7.81±0.44 e 2 fm 4 
13.4±1.8 e 2 fm 4 


8.8 e 2 fm 4 
19.3 e 2 fm 4 


8.0 e 2 fm 4 
3.5 e 2 fm 4 



the 3a cluster model, because (3 decay of nucleons in a clusters is forbidden. In other words, 
the origin of the (3 decay strength is the components of the dissociation of the a clusters. 
On the other hands, even in shell-model calculations with the large model space of (0 + 2)hu 
the O3" state does not appear in low energy region as mentioned above. Therefore there 
has been no theoretical approach which reproduce the strength of j3 decay to the excited 
+ states. This is the first microscopic calculation that can describe a large number of the 
experimental data of the transition strength to the excited states of 12 C. It means that in 
the present calculations the component of the dissociation of a clusters is reasonably large 
enough to explain the /3 transition strength. This success is due to the flexibility of AMD 
wave functions which can describe both shell-model-like and molecule-like aspects. 

Another shell- model-like aspect is found in the lf(T = 0) state which is considered to 
be a shell-model-like state and can not be described in 3a cluster models. In this work the 
(3 transition strength to if from 12 N reasonably agrees to the experimental data for the 
1 + (T = 0) state at 12.71 MeV, though the excitation energy is overestimated (see Table H 
and Fig [ID . 

In summary, we achieved the variation after the parity and the total angular momentum 
projection within the framework of antisymmetrized molecular dynamics and applied it to 
the excited states of 12 C. The framework is a new approach for the study of the structure of 
excited states with freedom from any model assumptions such as inert cores and the axial 
symmetry. It is a practical method of the variation after spin parity projection with finite 
range interactions. The model can represent both clustering structures and the dissociation 
of clusters due to shell-model-like aspects within one framework. This microscopic model is 
a powerful tool which is applicable to stable as well as unstable nuclei. Many excited levels 
of 12 C, E2 transition strengths and (3 decays to the excited states are reproduced within 
one microscopic framework for the first time. We can obtain information of the structures 
of many excited states from the experimental data of the transition strength. It owes to the 
flexibility of the AMD wave function which can describe both shell-model-like and clustering 
aspects seen in the excited states. 

The author would like to thank H. Horiuchi for many discussions and advises. She is also 
thankful to N. Itagaki, Y. Kondo and O. Morimatsu for helpful discussions. She is grateful to 
A. Ono for technical advices on numerical calculations. This work was partly supported by 
Research Center for Nuclear Physics (RCNP) of Osaka University and Institute of Physical 
and Chemical Research (RIKEN) for using computer facilities. 
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TABLE II. The experimental data for (3 decays 12 N(/3 + ) 12 C compared with the theoretical 
results. 



Decay to 12 C* (MeV) 


J n 


exp. 


theory (jf) 





0+ 


4.120 ± 0.003 


3.8 (0+) 


4.44 


2+ 


5.149 ± 0.007 


4.8 (2+) 


7.65 


0+ 


4.34 ± 0.06 


4.0 (0+) 


10.3 


(0+) 


4.36 ± 0.17 


4.7 (0+) 


12.71 


1+ 


3.52 ± 0.14 


3.8 (1+) 
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FIG. 1. Excitation energies of the levels of 12 C. In the right-hand side, VAP calculations in the 
AMD framework are shown and compared with the experimental data (left). 



6 



